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Abstract 


A  conraiunications  network  that  is  part  of  a  soft  real-time  system  will  need  to  transmit 
messages  within  a  bounded  delay,  but  may  allow  some  messages  to  miss  this  bound  and  be 
dropped  within  a  Tnaximiim  prespecified  rate  of  message  loss.  This  paper  presents  a  media  access 
protocol  for  soft  real-time  systems  implemented  on  a  slotted  radio  channel  with  binary  feedback. 
The  protocol  is  based  on  the  Gallager  FCFS  and  Capetanakis  splitting  algorithms,  but 
incorporates  strict  delay  bounds  using  packet  laxities.  Also  presented  is  an  analytic  model  for  this 
protocol  by  examining  the  probable  lengths  of  the  collision  resolution  intervals  given  the  current 
lag  at  any  time.  Both  analytic  and  simulation  results  are  obtained  to  study  the  maximum  input 
tr^c  rates  that  can  be  sustained  for  various  laxities,  delay  bounds,  and  message  loss  rates. 
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!•  Introduction 


Communication  on  the  battlefield  often  is  made  up  of  messages  with  deadlines.  If  a  message  is  not 
delivered  to  its  destination  on  time,  the  data  are  no  longer  useful.  Not  only  can  the  lost  information  have 
an  adverse  effect  on  battlefield  management,  but  current  communication  protocols  do  not  take  deadlines 
into  account,  resulting  in  transmission  of  out-of-date  messages  and  making  poor  use  of  the  channel.  Ideally, 
transmission  scheduling  will  take  into  account  both  the  needs  of  applications  as  well  as  individual  message 
deadlines.  On  the  battlefield,  where  bandwidth  is  limited,  this  is  especially  important. 

Battlefield  applications,  since  they  directly  interact  with  the  world,  often  must  respond  to  changes  within 
some  predetermined  amount  of  time.  Such  systems  are  classified  as  hard  or  soft  real-time  systems.  A  hard 
real-time  system  requires  that  all  deadlines  are  always  met  because  of  the  importance  of  the  data:  aeronautics 
or  nuclear  power  control,  for  example.  Soft  real-time  systems  can  afford  to  permit  some  deadlines  to  slip. 
The  battlefield  can  be  modeled  as  a  loosely  coupled  soft  real-time  system  because  each  node  functions 
independently  but  requires  high-level  coordination  between  nodes.  In  such  a  communication  network,  quality- 
of-service  constraints  may  require  that  a  message  be  transmitted  within  a  bounded  delay,  but  that  some 
messages  may  miss  this  bound  and  be  dropped,  resulting  in  a  certain  amount  of  message  loss.  For  instance, 
infrequently  dropped  voice  packets  do  not  interfere  with  a  conversation,  and  similarly,  tracking  objects  does 
not  require  processing  of  every  return  radar  signal.  In  either  the  hard  or  soft  case,  however,  all  layers  of  the 
protocol  stack  must  support  the  concept  of  deadlines.  Ultimately,  it  is  the  lowest  layer  that  provides  timely 
access  to  the  communications  medium.  In  this  paper,  a  media  access  protocol  for  soft  real-time  systems 
implemented  on  a  slotted  radio  channel  with  binary  feedback  is  considered. 

Because  each  application  has  its  own  requirements  for  a  minimum  acceptable  level  of  protocol  perfor¬ 
mance,  our  interest  is,  given  these  requirements,  to  determine  the  maximum  traffic  rate  the  system  can 
handle.  Two  application  requirements  are  considered:  minimum  fraction  of  packets  that  must  be  success¬ 
fully  transmitted  within  the  specified  delay  bound,  and  average  delay  experienced  by  a  packet.  Knowing 
these  for  various  values  of  system  parameters  will  allow  the  protocol  to  be  tailored  for  use  by  a  variety  of 
soft  real-time  applications  without  designing  to  the  lowest  common  denominator,  i.e.,  the  safest,  lowest  loss, 
and  lowest  throughput  system. 

The  obvious  disadvantage  to  using  a  general  random  access  scheme  for  real-time  communication  is  that 
the  worst  case  channel  access  time  is  unbounded  because  of  packet  collisions.  The  problem  is  further 
compounded  by  the  fact  that  general  algorithms  do  not  take  packet  transmission  deadlines  into  account. 

Approaches  to  limiting  or  removing  these  shortcomings  for  time-constrained  communication  have  concen¬ 
trated  on  two  methods:  the  use  of  virtual  time  clocks,  and  window-splitting  techniques.  Virtual  time  clocks 
were  first  proposed  by  Molle  and  Kleinrock  [1]  and  based  on  message  arrival  time.  This  has  the  advantage 
of  making  transmission  of  queued  messages  fairer.  The  method  was  adapted  by  Ramamritham  and  Zhao  [2] 
to  take  into  account  various  time-related  properties  of  a  packet  for  soft  real-time  systems  and  shown  via 
simulation  to  work  better  than  protocols  not  designed  for  real-time  use.  Subsequently,  Zhao,  Stankovic, 
and  Ramamritham  [3]  proposed  a  window-splitting  protocol  that  always  performed  in  simulation  at  least  as 
well  as  the  virtual  time  protocols  and  often  better.  Less  complex  window-splitting  algorithms  are  presented 
for  both  hard  and  soft  real-time  systems  by  Arvind  [4],  where  protocol  operations  are  simulated  and  some 
worst-case  performance  analysis  is  also  presented.  Paterakis,  Georgiadis,  and  Papantoni-Kazakos  [5]  present 
a  simple  protocol  appropriate  for  some  soft  real-time  systems  and  perform  an  in-depth  analysis  of  it.  Pater¬ 
akis’  technique  of  analysis  is  followed  in  this  report  and  extended  to  analyze  a  more  complicated  protocol  that 
would  also  be  suitable  for  use  in  a  soft  real-time  environment.  The  protocol,  the  CG  (Capetanakis-Gallager) 
protocol,  is  a  combination  of  the  Capetanakis  splitting  algorithm  and  the  Gallager  FCFS  algorithms  [6],  [7] 
with  binary  feedback  and  our  addition  of  strict  delay  bounds. 

2.  Algorithm 

In  this  protocol,  a  packet  has  a  single  property  of  interest,  its  laxity.  Laxity  is  the  maximum  amount 
of  time  that  can  elapse  prior  to  transmission,  after  which,  the  packet  will  not  reach  its  destination  on  time 
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(only  single-hop  radio  channels  are  considered  in  this  analysis).  Once  a  laxity  is  assigned  to  a  packet,  at  each 
tick  of  the  clock  it  is  decremented  and  the  packet  discarded  should  it  reach  zero.  The  channel  is  accessed 
in  a  slotted  manner  with  one-slot-long  packets  with  binary  (collision/noncollision)  feedback.  Collisions  are 
resolved  using  the  technique  described  below. 


2.1  CG  CRA 

When  two  or  more  nodes  transmit  at  once  and  collide,  the  collision  resolution  algorithm  (CRA)  com¬ 
mences  and  only  nodes  involved  in  the  collision  may  contend  for  the  channel.  Only  after  all  collided  packets 
have  been  either  successfully  transmitted  or  else  discarded  due  to  missed  deadlines,  can  other  nodes  again 
contend  for  the  channel.  During  the  collision  resolution  interval  (CRI),  a  window  of  initial  length  A  is  used, 
and  packets  within  it  are  ordered,  left  to  right,  from  lowest  laxity  to  highest.  In  the  slot  after  a  collision,  only 
packets  whose  laxities  fall  within  the  window  may  transmit.  If  another  collision  occurs,  the  window  is  split 
in  half,  and  the  CRA  recurses  first  on  the  left  half  and  then  on  the  right.  The  CRA  behaves  similarly  to  the 
classical  FCFS  splitting  algorithm  [6]  with  a  few  differences:  packets  in  a  window  following  a  collision  are 
ordered  by  laxity  rather  than  by  node  arrival  time,  and  feedback  is  binary  rather  than  ternary.  Additionally, 
each  window  is  treated  in  the  recursive  manner  of  the  Capetanakis  algorithm  rather  than  returning  right 
subwindows  to  the  waiting  interval  as  in  the  Gallager  algorithm.  That  is,  in  the  Capetanakis  algorithm,  at 
any  point  in  time,  there  can  be  an  arbitrary  number  of  windows  active.  In  the  Gallager  algorithm,  there 
are  never  more  than  two  windows  active  at  once.  A  further  constraint  is  a  bound  T  that  is  the  maximum 
number  of  slots  a  CRI  may  comprise.  If  T  slot  times  are  exceeded,  packets  involved  in  the  CRI  are  dropped, 
and  the  algorithm  resets. 

Because  of  the  laxity  ordering,  the  first  successful  transmission  during  a  CRI  will  be  the  lowest  laxity 
packet,  the  second  will  be  the  second  lowest,  and  so  on,  guaranteeing  a  laxity-ordered  transmission  schedule 
appropriate  in  a  real-time  setting.  During  a  CRI,  an  initial  collision  that  happened  some  while  ago  is  being 
resolved.  This  means  that  there  is  a  lag  of  d  units  between  “now”  and  the  period  of  time  currently  being 
examined  by  the  CRA.  The  lag  d  is  important  in  a  real-time  environment  because  the  lag  induced  by  the 
CRI  means  that  only  T-d  slots  remain  before  a  packet  with  initial  laxity  T  must  be  transmitted  or  dropped. 
In  addition,  the  width  u  of  a  window  plays  a  crucial  role  in  the  length  of  a  CRI.  Short  initial  windows  waste 
time  because  many  will  be  empty,  while  long  ones  potentially  encompass  several  packets  leading  to  still  more 
collisions. 

Letting  ft  denote  the  feedback  corresponding  to  slot  t:  /<  =  c  if  there  wcis  a  collision  in  slot  t,  otherwise 
ft  =  nc.  In  addition,  the  position  of  the  left  edge  of  the  window  at  time  t  is  denoted  as  xt  and  its  length,  in 
slots,  as  Uf.  Finally,  ht  can  take  on  the  value  of  L  or  indicating  whether  the  collision  resolution  algorithm 
was  in  the  left  or  right  sub  window.  By  convention,  the  algorithm  is  initially  in  the  right  half.  At  time  t  =  0 
the  system  is  empty  and  that  slot  has  a  feedback  value  of  /o  =  nc.  The  CG  CRA  operates  as  follows: 

1.  If  collision  resolution  is  not  in  progress,  then  a  node  with  a  packet  that  arrived  in  slot  t  —  1  may 
transmit  it  in  the  current  slot  f. 

2,  If  collision  resolution  is  in  progress: 

♦  If  ft-i  =  nc  and  ht^i  =  L,  then  =  Xt-i  +  at  =  Ut-i,  and  ht  =  R. 

♦  If  /t_i  =  nc  and  ht~i  =  J?,  then  xt  =  +  at_i,  =  min(2at_i,min(d.  A)),  and  ht  =  R. 

♦  If  ft-i  =  c,  then  Xt  =  at  =  and  ht  =  L. 

2.2  Paterakis  CRA 

Paterakis  et  al.  [5]  analyze  a  time-bounded  CRA  that  does  not  take  deadlines  into  account.  Their  protocol 
is  elegantly  simple.  The  protocol  is  described  here,  and  its  performance  is  later  compared  to  that  of  the  CG 
CRA.  Note  that  deadlines  are  not  taken  into  account  though  the  CRI  length  is  bounded.  Each  node  has  a 
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counter  whose  value  may  be  either  1  or  2,  but  can  transmit  only  when  the  counter  value  is  1.  Upon  collision, 
a  CRI  begins.  During  the  CRI,  after  each  collision  all  colliding  nodes  flip  coins.  Counters  are  assigned  a  new 
value  of  1  on,  say,  heads,  and  2  on  tails.  After  each  non  collision  slot,  all  nodes  in  the  CRI  set  the  counter  to 
1  and  transmit.  A  CRI  ends  when  there  are  two  consecutive  noncollision  slots.* 

A  node’s  counter  at  time  t  is  denoted  by  r^.  With  this  notation,  the  algorithm  can  be  more  rigorously 
described. 

1.  A  packet  is  successfully  transmitted  if  and  only  if  =  1  and  ft  =  nc. 

2,  The  transitions  are: 

•  If  ft^i  =  nc  and  =  2,  then  r*  =  1. 

•  If  ft-i  =  c  and  rt-i  =  2,  then  Vt  =  2. 

•  If  /t-i  =  c  and  =  1,  then  rt  =  1  with  probability  0.5  and  rt  =  2  with  probability  0.5. 


Figure  1.  Relation  between  some  variables  in  algorithm. 

Figure  1  illustrates  some  of  these  variables  and  how  they  are  related.  In  the  figure,  the  current  moment 
of  the  illustration  is  time  t  and  because  of  some  previous  collision,  the  current  lag  is  d  slots.  The  previous 
CRI  was  completed  at  time  ti.  As  indicated,  the  time  t2  corresponds  to  a  laxity  of  one  slot.  Anything  with  a 
smaller  laxity  cannot  be  transmitted  before  its  deadline  and  so  is  rejected,  or  discarded,  outright.  Similarly, 
to  successfully  transmit  all  packets  in  the  current  window,  the  CRI  must  complete  within  the  next  =  T  —  d 
slots.  The  variable  lu,d  will  be  defined  in  detail  in  subsequent  sections  but  is  simply  the  length  of  the  CRI 
that  begins  at  time  t  and  ends  at  time  t'. 


3,  Analysis 
3.1  Notation 

Our  analysis  technique  and  notation  are  based  on  the  methodology  presented  by  Paterakis  et  al.  [5]. 
As  mentioned,  two  considerations  vital  to  successful  transmission  of  a  packet  are  the  current  lag  d  and 
the  number  of  slots  u  to  be  examined.  Appropriately,  the  following  variables  are  subscripted  with  this 
information. 

Uu^d  •  Number  of  packets  in  window  of  size  u  with  lag  d  that  are 
successfully  transmitted  during  the  CRI. 

Zu^d  Sum  of  delays  after  time  t  of  packets  in  Uu^d- 

'ipu,d  •  Sum  of  delays  of  the  n^^d  packets  within  the  current  window  of  length  u. 
lu,d  •  Number  of  slots  needed  to  examine  u  slots  when  current  lag  is  d  slots. 

In  order  to  compute  the  expected  values  for  the  fraction  of  traffic  transmitted  within  the  laxity  bound  and 
for  the  delay  experienced  by  successful  transmissions,  the  aforementioned  variables  must  be  incorporated 

*Most  window-splitting  techniques  split  based  on  some  parsimeter  such  as  arrival  time,  laxity,  etc.  When  two  packets  He, 
i.e.,  have  the  same  value  for  that  parameter,  rcuidomness  is  introduced  to  artificially  separate  the  values.  Paterakis  simply  uses 
this  technique  at  the  outset. 
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into  expressions  that  reflect  not  just  what  occurs  for  given  window/lag  values  but  for  full  CRI  lengths.  The 
length  of  a  CRI  is  the  nunaber  of  slots  it  takes  to  move  from  a  lag  of  d  =  1  to  the  next  d  =  1  lag.  Therefore, 
the  subscript  d  on  the  following  variables  represents  what  happens  between  the  current  lag  of  d  slots  and 
when  the  lag  next  returns  to  1.  When  d  ^  1,  the  algorithm  is  at  some  intermediate  point  in  the  CRI. 

hd  :  Given  current  lag  of  d,  the  number  of  slots  needed  to  return  to  lag  of  1. 

Wd  ■  Cumulative  delay  of  packets  transmitted  during  the  hd  slots. 
ad  :  Number  of  packets  transmitted  during  the  hd  slots. 

Also,  given  that  some  event  occurs  during  a  CRI  of  length  /,  we  must  multiply  that  event  by  the  probability 
that  the  CRI  is  actually  of  that  length.  Thus,  the  form  of  the  final  type  of  variable  is  denoted  as: 

P(l  I  Ujd)  :  Given  that  u  slots  are  to  be  examined  and  the  lag  is  presently 
d  slots,  the  probability  that  the  CRI  will  be  /  slots  long. 

Finally,  to  compute  probabilities  of  a  given  number  of  packets  arriving  in  some  interval  u,  the  Poisson  arrival 
process  is  used. 


3.2  Model 

In  this  analysis,  a  simplifying  assumption  is  made,  namely  that  the  initial  laxities  of  all  packets  are 
identical  and  equal  to  T,  Because  of  this,  the  laxity  ordering  of  the  packets  is  in  fact  the  arrival  time  plus 
a  constant  laxity  offset.  Later  simulation  results  are  presented  for  the  case  where  initial  laxities  are  not 
identical. 

Considering  lags  of  1  slot  as  regeneration  points,  a  single  CRI  is  one  cycle  of  a  regenerative  stochastic 
process  as  pointed  out  by  Georgiadis,  Merakos,  and  Papantoni-Kazakos  [8].  Defining  Z  =  £*{^1}  as  the 
number  of  packets  successfully  transmitted  in  a  CRI,  and  H  =  £{^1}  as  the  expected  length  of  a  CRI,  then 
ZjH  is  the  traffic  rate  of  successful  packets.  This  rate  must  be  less  than  or  equal  to  the  original  rate  A. 
Therefore,  the  fraction  p  of  generated  packets  that  are  successfully  transmitted  is 

(1) 

If  we  next  define  W  =  £{u^i}  as  the  cumulative  expected  delay  of  all  packets  in  the  CRI,  then  the  per  packet 
expected  delay  is  simply 

0=f.  (2) 

In  the  following  sections,  values  for  H,  Z,  and  W  are  derived. 


3.3  Recursions 

In  the  following  sections,  let  us  denote  by  £{  •  |  fc,  £  }  the  expected  value,  given  that  h  packets  are  in  the 
current  window,  and  the  maximum  length  of  the  collision  resolution  interval  is  B.  Let  P(-  \k,B)  denote  the 
conditional  probability  with  given  variables  as  above. 

3.3.1  Probable  Length  of  CRI 

The  probability  of  a  CRI  lasting  /  slots  can  be  expressed  recursively.  Let  P{1  \  u,  d)  denote  the  probability 
that  a  CRI  is  I  slots  long  given  that  u  slots  must  be  examined,  and  that  the  current  lag  is  d  slots.  This 
probability  is  the  sum  of  the  probabilities  that  the  u  slots  contain  0,  1,  . . . ,  00  packets  and  that  these  packets 
can  be  successfully  transmitted  within  the  remaining  time,  T  -  d. 

=  (3) 

fc=0 
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Note  that  in  the  expression  above,  P{1  |  u^d)  is  in  terms  of  P{1  |  k^B)  where  k  is  the  number  of  packets, 
and  B  is  the  laxity  bound.  Initial  conditions  are  listed  below.  In  the  general  case  of  a  collision,  however,  the 
algorithm  states  that  the  current  window  will  be  split  in  two.  Therefore,  the  probability  that  a  CRI  is  of 
length  /  given  that  a  collision  occurs  and  wastes  one  slot,  is  the  probability  that  the  sum  of  the  CRI  lengths 
of  the  two  subwindows  is  /  —  1.  Of  the  k  packets  involved  in  the  collision,  the  probability  of  a  packet  being 
in  one  half  versus  the  other  is  simply  For  k  >2  and  S  >  2,  this  is  recursively  expressed  as 

P{l\k,B)  =  P(l^j\i^B^l)  (4) 

'  P{j  —  1  \  k  —  i,B  —  1  —  (I  —  j)),  with  probability  2“^ (^) . 

This  is  expressed,  for  /  >  1, 5  >  1,  as 

j=i i=0 

with  the  following  initial  conditions: 

P(0  |ifc,0)  =  1, 

P(0  I  jfe,B)  =  0,  B>0 

P(/|fc,0)  =  0,  />0 

P(/|jk,S)  =  0,  1>B 

P{l\k,B)  =  l,  B>l,0<k<l, 

P{l\k,l)  =  l,  k>2, 

P{l\k,B)  =  0,  B>2,k>2. 

Because  the  number  of  packets  k  and  the  laxity  T  are  finite,  the  probabilities  of  CRI  lengths  can  be  determined 
by  calculating  a  finite  number  of  terms. 


Q)  [Pii -j\i,B-l)-Pij-l\k-i,B-l-(l- j))] 


3.3.2  Expected  Length  of  CRI 

When  the  lag  d  is  such  that  A  <  d,  the  expected  length  of  a  CRI  is  determined  recursively  as  follows. 
With  a  lag  of  d,  there  are  d  slots  to  examine.  However,  only  A  slots  at  a  time  will  be  considered.  As  a  result, 
the  expected  CRI  length  ha  with  lag  d  is  the  number  of  slots  it  takes  to  examine  A  slots  plus  the  expected 
CRI  length  of  now  only  d  —  A  remaining  slots  plus  the  additional  number  of  slots  it  took  to  analyze  that 
first  window  of  length  A.  This,  as  well  as  when  the  lag  is  small  and  d  <  A,  is  expressed  as 

{  ld,d,  ld,d  =  1,  1  <  d  <  A 

hd=l  ld,d  +  hi^^^,  I  <  ld,d  <T  -  \d],  l<d<A  (6) 

[  /A,d  +  /ld-A+?A,<i>  A<d<T-l. 

The  expected  value  of  a  CRI  given  length  u  to  be  examined  and  lag  d  is  developed  shortly.  However,  any 
given  CRI  is  always  of  integer  length.  To  determine  the  expected  length  of  the  remainder  of  the  CRI  after  A 
slots  are  resolved  means  giving  consideration  to  each  of  the  possible  lengths  of  the  CRI  for  the  first  A  slots. 
The  results  of  Section  3.3.1  are  used  for  this,  keeping  in  mind  that  the  CRI  length  can  range  between  1  slot 
and  T  —  d  slots.  So,  taking  expectations  and  denoting  Hd  =  E{hd}  yields 

T-\d\ 

E{ld,d}+  E  ^mP(m|d,d),  l<d<A, 

H<i={  (7) 

E{lA,d}-\-  E  I  A, d),  A<d<T— 1. 

m=l 

This  system  of  equations  is  finite  because  of  the  bounded  nature  of  the  algorithm.  See  Paterakis,  Georgiadis 
and  Papantoni-Kazakos  [5]  for  further  mathematical  discussion.  The  system  is  represented  by  an  n  x  n-h  1 
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matrix  where  all  elements  are  the  conditional  probabilities  except  for  the  rightmost  column.  That  column 
is  made  up  of  (the  negatives  of)  the  expected  values. 

The  expected  values  of  CRI  lengths  are  derived  below.  Equation  6  above  is  taken  directly  from  Paterakis, 
Georgiadis,  and  Papantoni-Kazakos  [5]  since  it  is  a  general -description  of  any  bounded  CRA.  Equation  7, 
the  expected  value,  addresses  the  analytical  details  specific  to  the  bounded  CG  CRA.  So,  given  a  length  of 
u  time  units  to  resolve  with  a  current  lag  of  d  time  units,  the  expected  length  of  that  CRI  is  derived  by  first 
summing  against  all  possible  numbers  of  packet  arrivals  k  in  the  window  u  and  multiplying,  respectively,  by 
their  probabilities  of  occurrence  yielding 


E{lu,,)  =  ^E{h,i\k,T-\d\  ^ .  (8) 

fc=0 

Considering  each  term  of  the  summation  individually,  when  it  is  given  that  the  window  u  contains  k  packets 
and  there  is  a  laxity  bound  5,  initial  conditions  are  straightforward.  If  0  or  1  packets  are  in  the  window, 
the  CRI  length  will  be  one.  When  ib  >  2,  the  resulting  collision  takes  one  slot.  Furthermore,  though  the 
k  packets  are  each  equally  likely  to  be  in  either  the  left  or  right  subwindow,  the  sum  of  the  bounds  of  the 
two  windows  must  be  at  most  jB  —  1  since  the  collision  used  one  already.  The  left  subwindow  always  takes 
at  least  one  slot,  and  in  some  cases  can  even  use  all  5  —  1  slots.  Any  remaining  slots  of  the  original  5  —  1 
can  be  used  by  the  right  subwindow.  Because  a  given  CRI,  as  opposed  to  the  average  case,  can  be  only  an 
integer  number  of  slots  between  one  and  T  —  [cT] ,  each  possible  length  must  be  considered.  The  summation 
with  index  m  Equation  9  loops  through  the  range  of  possible  CRI  lengths,  recurses  on  each,  and  multiplies 
each  results  by  the  probability  as  given  in  Equation  5. 

Denoting  Lk,B  =  E{  lu,d  1^,5},  the  recursive  expression  for  the  length  of  a  CRI  based  on  the  operation 
of  the  algorithm  is 


Lk,B 


1,  0<k<l 

1  -1-  +  E  Ei'^  \hB-  l)Lk-i,B-i-m,  with  probability  k  >  1. 

m=l 


(9) 


For  initial  conditions  Lk^o  =  0  for  all  k^  and  Lq^b  —  =  1  for  5  >  1,  then 


Lk,B 


B-l 

1  +  Li^B^l  ^  ^  P{Tn  I  Z,  5  —  l)Lk-^i,B-l—m 

m—\ 


Note  that  the  bounded  nature  of  the  algorithm  guarantees  a  finite  expansion  of  terms. 


(10) 


3.3.3  Packets  Successfully  Transmitted 


The  next  value  to  be  determined  is  the  expected  number  of  successful  transmissions  between  times  with 
a  lag  of  1,  that  is,  a  single  CRI.  This  requires  knowledge  of  the  expected  length  of  a  CRI  as  developed  in  the 
preceding  section.  Given  a  current  lag  of  d  slots  with  d  >  A,  there  will  be  some  number  of  packets  transmitted 
during  examination  of  the  current  window.  This  examination  removes  A  slots  from  the  lag  d  but  introduces 
more  slots.  Thus,  the  total  expected  number  of  transmitted  packets  would  be 
also  account  for  the  situation  where  d  <  A,  meaning  the  lag  is  very  near  the  current  moment  in  time,  the 
following  expression  is  more  general: 


ad  = 


^djd  Idfd  —  Ij 

nd,d  +  0‘ii,i,  '^  <  U,d  <T  -  [d\, 

nA,d  +  ad-A+lA.d) 


Taking  expectations  and  denoting  Ad  =  E{a^, 


1  <  d<  A, 

1  <  d<  A, 
A<d<T-l. 


Ad={ 


T-\i\ 

E{nd^d}-k  E  AmP{m\d,d), 

m=2 

T-fdl 

m=l 


l<d<A, 

A<d<T-l. 


(11) 


(12) 
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.  Finally,  the  number  of  packets  transmitted  during  a  single  CRI  is 


E{nu.d}  =  f^E{  n„.d  \k,T-\d] 

A:=0 


(13) 


Denoting  Nk,B  =  E{nu,d  |  in  order  to  subsequently  calculate  the  expected  number  of  packets  suc¬ 

cessfully  transmitted  during  the  period  that  it  takes  to  move  from  a  lag  of  d  to  a  lag  of  1,  the  algorithm 
induces  the  relationship  below.  The  left  subwindow  can  use  up  to  JB  —  1  slots  to  transmit  the  i  packets  in  its 
window.  In  a  manner  identical  to  that  for  obtaining  the  expected  CRI  length,  each  possible  CRI  length  is 
considered  because  however  many  slots  a  given  left  window  CRI  requires,  the  right  window’s  CRI  has  that 
many  fewer  slots.  Because  individual  cases  are  being  considered,  the  expected  value  of  a  CRI  length  cannot 
be  used  as  a  subscript  on  the  last  factor  of  the  last  term. 


Nk,B 


0, 

1, 

J5-1 

+  E  \hB-  l)Nk-i^B-l-m, 


m—l 


.  For  initial  conditions,  =  0  for  all  and  Nq^b  =  0,  Ni^b  = 


2  =  0 


Ni,B-l  +  i 


m=l 


Ar=:0,B>  1, 
k  =  l,B>  1, 

k  >  IjWith  probability 
1  for  R  >  1,  then 

B  -l)Nk-i^B-l--^m  - 


(14) 


(15) 


As  before,  this  is  a  finite  expansion  because  of  the  properties  of  the  CRA. 


3.3.4  Delay  Experienced  by  Packets  Successfully  Transmitted 


Determining  the  cumulative  delay  experienced  by  successfully  transmitted  packets  is  perhaps  the  most 
complex  of  the  calculations  introduced  thus  far.  The  cumulative  delay  is  the  sum  of  the  delays  experienced 
within  the  window  u,  the  delays  experienced  by  the  current  lag,  and  the  delays  during  the  CRI.  That  is, 


-h  Zd4,  ld,d  =  1,  1  < 

+  Zd.^d  -h  ,  1  <  ld,d  <T  -  \d],  l<d<A, 

+  ZA,d 

^  +  (d- A)nA,d  +  u;^_A+?^.4,  A<d<r-1. 


(16) 


From  the  memoryless  nature  of  the  Poisson  process,  the  long-term  average  position  within  a  window  u  is 
the  middle.  Therefore  each  packet  in  the  window  experiences,  on  average,  a  delay  of  half  the  window  length. 
The  expected  value,  then,  is, 

Taking  expectations  where  Wd  =  E{wd}j  we  see  that 


Wd={ 


E{i^d,d  ^  Zd,d} WmB{m\d,d), 

m=2 

T-M 

E{'ipA,d  + ZA,d  + (d- A.)nA,d}  +  E  Wd-A+mP{m\  A,d), 

m=l 


1  <  d<  A, 
A<d<T-l. 


(17) 


.  Therefore,  the  cumulative  delay  of  packets  transmitted  during  a  single  CRI  is 

E{zu,d}  =  \k,T-\d\ 


(18) 


The  initial  conditions  are  easy:  a  window  with  no  packets  experiences  no  delay,  and  a  window  with  one  packet 
experiences  a  delay  of  one.  It  is  more  involved,  however,  when  a  collision  has  occurred.  After  splitting,  each 
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packet  in  the  left  window  experiences  an  additional  delay  of  one  slot  due  to  the  collision.  The  packets  in  the 
right  window,  though,  will  experience  the  total  delay  of  the  left  window’s  resolution  in  addition  to  a  similar 
one-slot  delay/packet  for  the  initial  collision  that  caused  the  split.  Similar  to  the  development  of  previous 
expressions,  the  number  of  slots  used  by  the  left  subwindow  is  subtracted  from  the  maximum  available  slots 
for  use  by  the  right  subwindow.  This  is  expressed  mathematically  as 


r  0, 

1, 

^  P{m  \  i,B  —  1) 

m—1 

^  [(m  +  l)Nii  +  — mjj 

With  the  initial  conditions  above,  i.e.,  =  0,  Zq^b  =  0,  and  Zi^b 


k  =  0, 
k=l, 

k>2,  with 
probability 
=  1. 


Zk^B  =  ^  ^Ni^B-1  +  Zi^B-1  +  ^  Pim  \i,B-  l)[{m+l)Nk-i,B-l-m  +  Zk-i,B-\-m^  ■ 


(19) 


(20) 


This,  again,  is  a  finite  expansion. 


4.  Evaluation 

In  real-time  systems  it  is  important  that  the  application’s  requirements  be  met  by  all  layers  of  the  protocol 
stack.  As  previously  mentioned,  at  the  media  access  layer  of  a  soft  real-time  system,  this  translates  to  a 
guarantee  that  some  minimum  fraction  of  traffic  is  in  fact  transmitted.  Similarly,  while  there  is  a  maximum 
laxity  value,  it  is  also  desirable  to  design  to  a  significantly  smaller  average  delay.  Continuing  to  follow  the 
notation  of  Paterakis  et  al.  [5],  these  design  parameters  are  called  ei,  for  the  fraction  of  packets  transmitted 
successfully,  and  e2,  for  the  delay  experienced  by  an  average  packet. 

While  it  is  the  application  that  drives  the  choices  for  maximum  laxity  T  and  for  ei  and  62,  the  initial 
window  width  A  is  a  basic  design  parameter  of  the  protocol  itself.  That  is,  it  is  chosen  from  the  optimization 
of  the  analysis  of  the  previous  section.  The  window  can  range  between  zero  and  T  —  1.  For  this  range  we 
would  like  to  find  the  maximum  system  traffic  rate  A*  that  allows  successful  transmission  of  at  least  ei  of 
the  traffic.  Because  of  the  complexity  of  the  expressions  developed,  analytic  optimization  is  difficult.  For 
given  values  of  A,  however,  the  problem  is  simple  to  solve  numerically  and  reduces  to: 

Ay^^,  =sup(A  :  /?(r(A,A)  >  ei).  (21) 

Because  of  space  considerations,  results  are  presented  only  for  A  =  3,  though  data  have  been  generated  for 
various  A  values  that,  through  observation,  encompass  the  maximum  throughputs. 

A  simulation  study  of  this  protocol  was  also  conducted  using  Opnet,  a  network  simulation  package.  The 
simulations  accurately  modeled  the  overall  system  functionality  including  the  channel  behavior,  but  did  not 
explicitly  model  individual  stations.  Each  simulation  was  run  with  95%  confidence  intervals  that  the  fraction 
of  successful  transmissions  p  wets  within  ±0.005  of  the  steady-state  value.  Input  traffic  rates,  in  packets/slot, 
ranged  from  0.050  to  0.600  in  increments  of  0.005.  When  window  size  A  =  3.0  and  laxity  T  =  20,  Figures  2 
and  3  show  overlaid  graphs  of  analytical  and  simulated  results.  Because  of  the  close  correlation,  subsequent 
graphs  do  not  include  simulated  data  for  easier  readability.  Simulation  results,  however,  have  been  obtained 
for  all  presented  graphs  and  show  similar  close  correlations. 

In  Figure  4,  three  curves  are  graphed  showing  the  maximum  sustainable  input  traffic  rates  attainable 
when  window  size  A  =  3.0  and  the  indicated  values  of  ei  are  used.  It  can  be  seen  from  these  figures  that 
the  maximum  input  rate  increases  as  initial  laxity  is  increased  and  also  as  the  threshold  for  success  rate  is 
decreased. 

Given  a  minimum  success  fraction,  an  additional  constraint  can  be  added  to  the  problem  posed  in 
Equation  21  so  that  an  average  delay  is  met  in  addition  to  meeting  the  maximum  laxity  T,  yielding: 

61,62  =sup(A  :  pt(A,A)  >  ei,J3T(A,A)  <  62).  (22) 
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Figure  2.  Simulated  and  analytical  results  for  fraction  successful  when  A  =  3.0  and  T  =  20. 
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Figure  3.  Simulated  and  analytical  results  for  average  delay  when  A  =  3.0  and  T  =  20. 
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Figure  4.  Maximum  input  traffic  rate  vs.  initial  laxity  with  minimum  success  constraints. 

When  an  average  delay  of  62  =  3  slots  is  met,  the  results  are  shown  in  Figures  5.  Similar  graphs  are  shown 
for  5  and  7  slots  in  Figures  6  and  7,  respectively.  All  three  figures  represent  performance  whh  window 
A  =  3.0  and  use  success  rate  constraints  of  ei  =  0.99,0.95,0.90.  It  is  interesting  to  note  that  with  the  new 
constraint  62)  the  curves  on  each  of  these  graphs  converge  whereupon  the  function  slowly  decreases.  This  is 
expected  when  increasing  traffic  rate  but  not  relaxing  the  delay  constraint. 


Figure  5.  Maximum  input  traffic  rate  vs.  initial  leixity  using  ei  and  62  —  3. 

It  could  be  expected  that  because  the  CG  CRA  makes  decisions  based  on  packet  deadlines  that  it 
would  outperform  the  Paterakis  CRA.  However,  the  CG  CRA  has  several  shortcomings  when  used  in  an 
environment  where  packets  have  fixed  initial  laxities.  Packets  with  very  close  laxity  values  cannot  be  quickly 
separated,  so  that  the  CRA  must  continually  split  until  the  windows  are  small  enough  to  separate  the 
packets.  Not  only  is  this  time  consuming,  but  the  CRA  must  then  recursively  double  the  window  lengths 
until  again  reaching  the  original  length  A.  Because  of  the  time-bounded  nature  of  the  CRA,  the  CRI  often 
expires  before  packets  with  similar  laxity  values  can  be  transmitted.  On  the  other  hand,  the  coin  flipping 
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Figure  6.  Maximum  input  traffic  rate  vs.  initial  laxity  using  ei  and  ej 


Figure  7.  Maximum  input  traffic  rate  vs.  initial  laxity  using  ei  and  62 


of  the  Paterakis  CRA  introduces  randomness  at  each  step  of  a  CRI,  making  it  more  likely  that  packets  will 
be  more  quickly  separated  into  two  groups.  Figure  8  illustrates  the  performance  differences.  It  shows  two 
sets  of  three  curves;  one  set  is  the  performance  of  the  Paterakis  and  the  other  of  the  CG  CRA,  representing 
maximum  sustainable  offered  loads  for  ei  =  0.99, 0.95, 0.90.  The  Paterakis  CRA  always  outperforms  the  CG 
CRA,  even  at  low  traffic  rates. 


Figure  8.  Maximum  input  traffic  rates  vs.  initial  laxity  with  minimum  success  constraints  using  only  ei* 


There  are  often  cases  where  soft  real-time  systems  generate  packets  with  variability  in  initial  laxities  unlike 
the  simply ing  assumption  made  in  this  analysis.  Using  the  simulation  model  to  study  such  a  situation,  each 
packet  was  assigned  an  initial  laxity  uniformly  distributed  in  the  range  [2,  T].  Figure  9  presents  these  results, 
showing  the  maximum  input  rate  achieved  as  the  average  laxity  is  increased.  Because  of  the  variability ,  the 
maximum  input  rates  are  now  smaller  than  the  rates  achieved  for  constant  laxities.  Figure  10  shows  the 
average  delays  experienced  for  the  corresponding  success  rates.  The  results  shown  represent  performance  for 
a  window  size  of  3.0  slots.  There  are  no  observable  differences  when  a  constraint  of  €2  is  imposed. 


Figure  9.  Maximum  input  traffic  rate  vs.  average  laxity. 


Average  Laxity  (slots) 


Figure  10.  Maximum  input  traffic  rate  vs.  average  laxity. 


5.  Conclusions 


To  support  a  soft  real-time  system,  a  network  must  incorporate  the  concept  of  deadlines  in  all  layers  of 
the  protocol  stack.  This  has  been  traditionally  difficult  to  do  in  random  access  protocols,  where  packet  delays 
may  be  potentially  unbounded  due  to  collisions.  By  assuming  that  a  message  is  dropped  whenever  its  delay 
exceeds  its  specified  laxity,  and  by  using  a  minimum  rate  of  successful  transmission,  the  analysis  presented 
has  shown  that  window-splitting  protocols  can  be  modified  to  work  successfully  in  these  environments. 
An  analytic  model  of  such  a  protocol  has  been  presented  that  can  be  used  to  determine  proper  operating 
parameters  for  specified  quality-of-service  constraints. 

A  comparison  of  the  CG  CRA  with  the  Paterakis  CRA  yielded  surprising  results.  It  showed  that  even 
though  the  Paterakis  CRA  does  not  take  deadlines  into  account,  its  extreme  simplicity  allows  it  to  perform 
better  than  the  CG  CRA.  It  seems  likely  that  a  variant  of  the  Gallager  CRA  would  perform  similarly  to  the 
Paterakis  CRA.  It  would  also  be  interesting  to  see  if  a  hybrid  protocol  can  be  developed  including  features 
of  both  the  Gallager  and  Paterakis  CRAs.  Such  a  protocol  could  be  of  significant  benefit  to  battlefield 
communication.  The  protocol  could  be  tailored  to  the  needs  of  applications  while  still  basing  transmissions 
on  message  deadlines. 
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